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Chapter  I 


Introduction 

J.R.  Heirtzler  and  E.G.  Zurflueh 

In  a  sense  this  is  a  "housekeeping”  report  in  that 
it  brings  together,  in  a  more  or  less  orderly  fashion,  for  our 
own  use  some  of  the  analysis  techniques  that  have  been  used  by 
various  persons  at  the  Lamont  Geological  Observatory.  It  is 
felt  that  such  material  may  be  of  interest  to  individuals  in¬ 
itiating  an  analysis  of  the  geomagnetic  field  intensity  for 
purposes  of  determining  earth  structure. 

In  the  measurements  of  the  geomagnetic  field  intensity 
in  marine  areas  by  Lamont  it  has  seldom  been  possible  to  take 
measurements  over  a  systematic  grid  and  construct  magnetic  contour 
charts.  Measurements  usually  consist  of  a  magnetic  profile 
along  the  track  of  a  research  vessel,  these  measurements  being 
taken  with  a  towed  magnetometer  that  records  the  total  intensity 
of  the  geomagnetic  field.  Over  a  period  of  time  our  ships  and 
other  ships  and  aircraft  have  made  sufficient  measurements  in 
some  areas  to  permit  a  general  description  of  the  magnetic  a- 

r 

nomalies  even  though  track  spacing  and  heading  prohibit  contour- 

i 

ing.  Frequently  such  a  general  description  can  substantiate  or 
refute  geophysical  hypotheses  and  suggest  guide-lines  for  future 
investigations. 

An  attempt  to  deduce  exact  geophysical  information 
from  a  limited  magnetic  profile  is  recognized  as  a  highly  un¬ 
certain  undertaking,  especially  when  control  of  the  time 


1-1 


1-2 


variations  is  doubtful.  Alternatively,  analyzing  a  high  densi¬ 
ty  of  profiles  on  a  statistical  basis  may  be  of  value  if  no 
errors  in  basic  physics  are  committed  and  the  limitations  of 
the  techniques  are  understood. 

The  structures  causing  magnetic  anomalies  can  fre¬ 
quently  be  treated  as  cylinderical,  or  two-dimensional,  i.e. 
extending  to  plus  and  minus  infinity  in  a  direction  parallel 
to  a  coordinate  axis.  Two-dimensional  structure  seems  to  be 
particularly  common  in  certain  marine  areas,  see  for  example 
(Vacquier,  et.al.,  1961)  and  (Heirtzler,  et.al.,  1962).  It  is 
instructive  in  such  cases  to  determine  some  of  the  configu- 
ations  of  magnetic  materials  that  can  produce  anomalies  similar 
to  those  observed.  An  electronic  digital  computer  program  to 
calculate  the  anomaly  that  would  be  caused  by  an  assumed  two- 
dimensional  structure  has  been  used  by  this  institution  for 
several  years.  It  was  originally  written  as  a  program  for  the 
IBM  650  by  M.  Landlsman  of  Lamont,  It  was  written  again  for 
the  IBM  1620,  because  of  the  convenience  of  that  computer.  It 
Is  the  formulation  of  this  second  program  that  Is  given  here 
but,  naturally,  both  programs  give  the  same  answers. 

In  order  to  give  some  crude  feeling  of  how  the  various 
geometrical  factors  affect  the  total  field  anomalies  a  brief 
selection  of  anomalies  for  simple  bodies  is  given  in  Chapter  III. 
With  the  large  number  of  geometrical  parameters  that  are  important 
in  this  calculation  It  is  clearly  impossible  to  give  a  complete 
handbook  of  anomalies  with  one  parameter  varying  at  a  time.  The 
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relative  importance  of  the  parameters  is  discussed  in  Chapters 
III  and  IV. 

In  Chapter  IV  a  close  examination  is  made  of  the  calcu¬ 
lated  anomalies  for  those  bodies  with  vertical  sides.  A  working 
set  of  rules  is  presented  to  determine  the  configuration  of  such 
bodies  from  observed  profiles.  In  that  Chapter  it  is  also 
shown  how  far  mathematical  theory  can  be  carried  without  having 
to  make  very  specialized  assumptions  about  body  configuration. 

In  connection  with  the  presentation  of  lnterpretational 
techniques  in  this  paper,  a  brief  listing  of  some  of  the  pertinent 
literature  is  given  here  with  short  comments  on  the  contents 
of  each  publication. 

Several  authors  give  sets  of  model  anomalies,  usually- 
accompanied  by  rules  for  depth  determination  and  formulas  for 
the  particular  model  used.  In  this  way  Haalck  (1927)*  Gulatee 
(1938)  and  Holland  (1948)  compute  model  anomalies  of  the  hori¬ 
zontal  and  vertical  components  for  different  bodies.  Nettleton 
(1942)  determines  the  vertical  intensity  anomaly  for  various 
geometric  configurations.  Henderson  and  Zietz  (1948)  treat  the 
effect  of  point  and  line  sources  with  respect  to  total  intensity 
measurements.  In  19£8  the  same  authors  published  a  paper  on 
the  use  of  magnetic  doublets  in  the  interpretation  of  total  in¬ 
tensity  anomalies.  Vacquier,  et.al.  (1951)  interpret  total  in¬ 
tensity  maps  with  the  help  of  three-dimensional  (prismatic) 
models  and  computed  second  vertical  derivatives.  Zietz  and 
Henderson  (1954)  use  three-dimensional  double  layer  model  fields 
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for  the  analysis  of  total  intensity  data.  Baranov  (1955) 
standardizes  his  model  calculations  by  reducing  measurements 
for  different  inclinations  to  vertical  field  conditions. 

Smellie  (1956)  gives  total  intensity  depth  factors  for  point 
and  line  configurations  of  poles  and  dipoles.  Prom  general 
mathematical  considerations  Smith  (1959  and  1961)  derives  in¬ 
equalities  which  can  be  used  to  estimate  depth  and  intensity 
of  magnetization  for  arbitary  bodies. 

Before  the  advent  of  the  airborne  magnetometer  the 
horizontal  or  vertical  components  were  measured  almost  ex¬ 
clusively  and  the  model  anomalies  were  treated  accordingly. 

Since  about  1946  total  intensity  measurements  are  used  in 
general  and  the  newer  literature  usually  refers  to  total  in¬ 
tensity  anomalies.  According  to  potential  theory  it  is  possible 
to  calculate  any  component  from  measurements  of  any  other  com¬ 
ponent  if  certain  conditions  are  fulfilled  (for  example  see 
symmetry  relations  on  page  2-9).  A  group  of  papers  deal  with 
the  conversion  of  one  component  into  another  and  with  the  re¬ 
lated  problems  of  computing  derivatives  and  upward  and  downward 
continuation. 

Skeels  (1947)  was  the  first  to  explicitly  point  out 
the  general  possibility  of  the  above  mentioned  conversions  by 
surface  integration.  Vestine  (1941)  and  Vestine  and  Davids 
(1945)  also  point  out  the  importance  of  integral  methods  in  mag¬ 
netic  interpretation.  Hughes  and  Pondrora  (1947)  describe  the 
calculation  of  the  vertical  intensity  out  of  total  intensity 
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data.  Skaels  and  Watson  (1949)  treat  the  problem  of  conversions 
In  a  more  general  way.  Henderson  and  Zletz  (1949*0  give  a 
method  for  the  computation  of  second  vertical  derivatives  and 
in  another  paper  (1949b)  they  analyze  the  upward  continuation 
of  the  field.  Peters  (1949)  treats  the  problems  of  upward  and 
downward  continuation  and  the  calculation  of  derivatives.  He 
also  derives  methods  for  the  direct  determination  of  basement 
structure.  Henderson  (i960)  relates  the  different  operations 
of  conversion  to  one  basic  formula  and  uses  expressions  suita¬ 
ble  for  use  on  electronic  computers. 

Another  group  of  papers  is  concerned  with  the  inverse 
problem  of  the  calculation  of  magnetic  anomalies  caused  by 
bodies  of  any  given  shape.  Gassmann  (1951)  describes  a  graphi¬ 
cal  method  of  effecting  the  triple  integration  necessary  for 
the  solution  of  this  problem.  Baranov  (1953)  works  out  another 
graphical  method  with  a  different  sequence  of  Integrations. 
Henderson  and  Zletz  (1957)  simplify  the  computation  of  a  total 
intensity  anomaly  by  projecting  a  map  of  a  given  body  in  the 
direction  of  the  earth's  field  and  using  a  modified  Gassmann 
process  for  the  integration.  Talwani  and  Ewing  (i960)  calculate 
the  anomaly  by  dividing  the  body  into  parallel  and  polygonal 
laminae*  similar  to  the  treatment  in  this  report.  Their  method 
is  put  in  a  form  which  can  be  used  for  digital  computers. 

Vajk  (1951)  and  Nettleton  (1954)  analyze  the  first 
step  in  magnetic  anomaly  interpretation,  namely  the  separation 
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of  anomalies  from  an  assumed  regional  field.  A  summary  of  the 
different  interpretation  methods  is  given  by  Kohler  (1958). 

His  work  also  contains  a  rather  complete  list  of  references 
which  includes  papers  from  the  large  Russian  literature  on  this 
subject.  Lastly,  there  are  extensive,  company  classified  tech¬ 
niques  used  by  the  geophysical  exploration  companies. 


Chapter  II 


The  mathematical  Expression  for  the  Magnetic  Anomaly 
over  a  Two-Dimensional  Body  of  Polygonal  Cross  Section 

U.  Talwani  and  J.R.  Heirtzler 

A.  Derivation  of  the  Basic  Formulas 

It  is  characteristic  of  geophysics  texts  to  derive  the 
expressions  for  the  horizontal  and  vertical  magnetic  anomalies 
by  analogy  with  the  expressions  for  gravitational  anomalies.  In 
the  present  report  the  expressions  for  the  horizontal,  vertical, 
and  total  field  anomalies  are  derived  from  basic  magnetostatic 
theory. 

Consider  the  volume  element  located  in 

an  external  magnetic  field  (Pig.  2-1). 
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We  assume  that  the  volume  element  has  a  uniform  mag¬ 
netization  M.  Let  ju>  be  the  magnetic  moment  of  the  element* 
Then 


The  magnetic  potential*  at  the  origin*  due  to  the  volume 

element  is 


P  - 


/•*„*  •*-  a 

*\)#/* 


AY-  At 


where 

fP  =  T*  *-j  V  *#■  *  *  . 


At  the  origin  the  potential  of  a  rod  of  cross  section 
Ax&t  and  infinite  in  the  y-direction  is 


r  «  c.x.  a*  f 

(X*-  *■  yL-*  2  V  ,/u  y 


l  *x  / 


The  vertical  magnetic  field  strength,  V,  is 


K  =  -  If  -  A*  /  ***"«-  *!<*'■*') 

3*  l  + 


) 


(2-1) 
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and  the  horizontal  field  strength,  H,  measured  in  the  x  direction 
is 


(  ~  /  (2-2) 

The  horizontal  field  strength  measured  in  the  y-direction  is  zero 
since  y  is  not  contained  in  the  expression  for  the  potential.  If 
M  is  the  induced  magnetization  this  means,  physically,  that  we 
have  not  included  demagnetization  effects.  The  above  represent 
the  magnetic  anomalies  caused  by  the  rod,  that  is  the  regional 
field  strength  is  not  included. 

Consider  the  lamina  shown  in  Pig.  2-2,  infinite  In  the 
y-direction  and  infinite  In  the  positive  x-direction. 


-i wszzm/:. 


At 


¥ 

*■ 

Pig.  2-2 
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By  the  use  of  eqs  (2-1)  and  (2-2)  one  gets 


~  (X*  <- 


i  ^a*  i 

t  X1--***  ) 


(2-3) 


H  *=  24*  J 


~0^-^  )u 


*  2. 
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+-  H-e  ' 

X',4-  *v 


(2-4) 


The  variables  x  and  z  in  (2-3)  and  (2-4)  represent  the  coordinates 
of  the  end  of  the  lamina.  If  we  had  taken  the  lamina  running  in 
the  opposite  direction,  as  indicated  in  Pig.  2-3 


Pig.  2-3 


wo  would,  have  obtained 


V  = 


4 

t  { 

^  K*--4-  ) 


H  =  *6* 


|  ~M«*  ~***L 
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Notice  that  V  and  H  bear  the  opposite  signs  to  eqs  (2-3)  an< 
(2— if) •  Next  consider  the  prism  shown  in  Fig*  2-4*  infinite 
the  y  direction* 


0 


i 
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the  integration  with  respect  to  z  in  eqs  (2-3)  and  (2-4): 


V 


H**-  H,* 

-i- 


A-t 


(2-5) 


+4 


Mx )i  •+*  t“1 

+  4X 


ib 


(2-6) 


The  integration  is  carried  out  from  the  small  to  the  larger  ©  . 
In  eqa  (2-3)  and  (2-4)  the  variables  x  and  z  represented  the  co¬ 
ordinates  of  the  end  of  the  lamina.  For  the  prism  x  and  z  are 
related  by  the  equation  for  the  line  AD: 


X  a  Cx,  f  l,  vck  $  )  - 


If  eqs  (2-5)  and  (2-6)  aw  integrated  and  the  variables  r  and 
used  instead  of  x  and  z  one  obtains: 


4 
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V  *  2.  tJ*  y  s*i-  p-Loj(rt  /r*,)| 

-  ^<^^(r*./r,)| j  (2  ^ 


H  -  X  Juu.  ^ 4>-c4>s<t>  ‘Cofiv't  /r,  )  J 
~  f+su.p  -Uf  (n/r,  J  jj \J 

(2-8) 

Had  it  bean  assumed  that  the  prism  was  oriented  as  in 
Pig.  2-5  the  above  would  have  been  preceded  by  a  minus  sign. 


0 


I 


2-9 


In  Pigs.  (2-i(.)  and.  (2-5)  the  prism  extended  infinitely 
in  one  of  the  x  directions.  It  is  clear  that  H  and  V  for  a  body 
of  finite  x  dimensions  can  be  treated  by  assuming  that  it  consists 
of  two  bodies  (ref.  Pig.  2-6).  One  must  be  careful  with  directions 
of  angular  measurements  and  with  signs  in  working  with  these 
equations. 


Pig.  2-6 


The  procedure  indicated  above  can  be  extended  to  any 
number  of  prismatic  or  polygonal  bodies  each  with  many  faces  and 
each  with  its  own  magnetization.  This  derivation  does  not  take 
into  account  nonuniform  magnetization  such  as  would  exist  near 
the  corners  of  the  bodies.  In  fact  the  actual  geological  bodies 
probably  do  not  have  distinct  comers  and  nonuniform  magneti¬ 
zation  is  believed  to  be  unimportant. 

It  is  interesting  to  note  the  symmetrical  relationships 
of  the  factors  in  eqs  (2-7)  and  (2-8).  That  symmetry  shows  that 
vertical  magnetization  effects  the  horizontal  anomaly  to  the  same 
extent  that  the  transverse  magnetization  effects  the  vertical 
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anomaly  and  that  the  transverse  magnetization  effects  the  hori¬ 
zontal  anomaly  to  the  same  extent  that  the  vertioal  magnetization 
effects  the  vertical  anomaly. 

B»  Induced  and  Remanent  Magnetization 
and  the  Total  Intensity  Anomaly 

If  one  wishes  to  utilize  eqs  (2-7)  and  (2-8)  for  the 

calculation  of  the  V  and  H  anomalies  due  to  induced  magnetization 

only  he  would  write 

M  =  k  F 

where  k  **  magnetic  susceptibility 

F  *  undisturbed  regional  total  intensity  vector 
and  as  illustrated  in  Pig.  2-7 
Px  =  k  P  cos  I  sin  s 
Mj  =  k  Pj  B  k  P  sin  I 


Pig.  2-7 


2-11 


I  «  inclination,  positive  if  F  is  belov  horizontal 
(northern  hemisphere) 

s  e  angle  of  strike,  measured  from  horizontal  pro¬ 
jection  of  F  and  either  positive  or  negative  y- 
axis.  s  must  always  be  less  than  180°  (and  thus 
sin  s  positive)  since  fill  specify  that  the  co¬ 
ordinate  axes  are  oriented  to  make  Fj  positive. 
This  orientation  of  the  axes  must  be  borne  in  mind 
for  the  oorrect  interpretation  of  results. 

If  one  wishes  to  utilize  eqs  (2-7)  and  (2-8)  for  the 
calculation  of  the  V  and  E  anomalies  due  to  a  body  with  remanent 
magnetization,  and  not  to  include  induced  anomalies,  he  would 
write 

Mx  *  Mrem  cos  a  sin  b 

M*  =  l*rem  si&  a 

where 

a  «=  inclination  of  remanent  magnetization  vector  below 
horizontal 

b  =  strike  angle  of  remanent  magnetization  vector 
measured  between  horizontal  projection  of  Mj,em 
and  positive  or  negative  y  axis,  less  than  180*. 

To  apply  the  equations  as  written  the  axes  must  be 
oriented  to  make  Mx  positive. 

If  it  is  desired  to  use  eqs  (2-7)  and  (2-8)  to  get  the 
V  and  H  anomalies  due  to  the  total  magnetization  (induced  and 
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remanent)  one  should  write 

“x  =  “tot  008  84  si"  p 
My  =  Mtot  sin  « 

where 

oc  =  inclination  of  total  magnetization  vector  below  hori¬ 
zontal 

fi  =  strike  angle  of  total  magnetization  vector  measured 
between  horizontal  projection  of  M^0t  and  positive  or 
negative  y  axis,  less  than  180*  To  apply  the  equations 
as  written,  the  axes  must  be  oriented  to  make  Mg  positive. 
Regardless  of  the  assumptions  used  to  determine  V  and 
H,  the  anomaly  in  total  intensity  is  obtained  as  follows: 

Let  T  be  the  anomaly  in  total  Intensity,  then 


(P  +  T)2  s  (Fx  +  H)2  +  (Py)2  +  (Fz  +  v)2 


Since  there  is  no  anomaly  in  the  y  direction  (ref.  p.  2-3 )• 
Squaring,  dropping  terms  that  contain  the  square  of  V  or  H  and 
remembering  that 

Fl  _ i.  i  i 

-  F. 


get 


Ref e ring  to  Pig.  2-7 


T  =  H  cos  I  sin  s  +  V  sin  I 


(2-9) 
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An  electronic  digital  oomputer  program  to  calculate 
V#  H,  and  T  is  given  in  the  Appendix. 

Eq  (2-9)  illustrates  that,  if  the  anomalies  are  small 
relative  to  F,  then  the  anomaly  in  total  intensity  is  the  com¬ 
ponent  of  V  and  H  in  the  Direction  of  F.  See  Fig.  2-8. 


Fig.  2-8 

The  magnitude  of  T  is  composed  of  the  projection  of  H 
(H  cos  I  sin  s)  and  of  V  (V  sin  I)  in  the  direction  of  F.  This 
geometrical  picture  is  invalid  if  the  anomaly  is  large  (we 
dropped  squared  terms  in  the  derivation  of  eq  2-9)  and  a  more 
adequate  vector  diagram  must  be  used. 
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Computed  Magnetic  Total  Intensity  Curves 
over  Two-Dimensional  Bodies 

G.  Peter 

A.  Description  of  the  Model 

With  the  aid  of  the  computer  program  (see  Appendix)  a 
few  families  of  magnetic  total  intensity  curves  have  been  com¬ 
puted  in  order  to  show  the  changes  in  these  curves,  when  the 
ambient  magnetic  field  directions  and  the  geometric  properties 
of  the  model  are  varied.  The  different  sets  of  curves  have  been 
obtained  by  varying  one  parameter  at  a  time.  Grouping  the  curves 
in  this  manner  enables  one  to  make  a  qualitative  analysis  of  them. 
The  curves  are  discussed  in  section  B. 


Pig.  3-1:  Ground  Plan  and  vertical 
section  of  assumed  two-dimensional  model  < 

x  =  line  of  measurement  s  =  strike 

d  =  depth  of  burial  W  =  Width 

D  =  vertical  extent 


The  model  used  in  these  calculations  is  a  two-dimension¬ 
al  rectangular  body,  which  has  infinite  length  in  a  horizontal 
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direction  (see  fig.  3-1) •  It  has  vertical  sides,  except  for  the 
group  6  curves,  where  the  sides  have  ±  30  and  ±  60  degree  angles 
from  the  vertical.  For  all  calculations  the  model  has  uniform 
isotropic  induced  magnetization  of  0.001  emu,  it  is  in  the  northern 
hemisphere  where  the  field  strength  is  50000  gammas.  The  pro¬ 
jection  of  the  magnetic  north  on  the  profiles  points  toward  the 
right . 


The  total  intensity  values  are  computed  for  a 


measurement  (profile),  which 
the  axis  of  the  body. 

The  various  curves 

1.  Variable  inclination 

2.  Variable  strike 

3.  Variable  depth 
Variable  width 

5.  Variable  vertical  extent 

6.  Slanted  bodies 


is  horizontal  and  at  right 

are  grouped  as  follows: 

I  (figure  3-2) 
s  (figures  3-3  to  3-7) 
d  (figures  3-8  to  3-12) 

W  (figures  3-13  to  3-15) 
D  (figures  3-1&  to  3-18) 
(figures  3-19  to  3-22) 


line  of 
angle  to 


These  curves  are  only  a  selection  for  a  few  varied 
parameters.  Additional  information  can  be  found  in  the  papers 
of  Haalck  (1927),  Gulatee  (1938)  and  Heiland  (1946). 


B.  Discussion  of  the  Computed  Curves 

The  following  is  a  short  discussion  of  the  general 
characteristics  of  the  computed  curves  and  of  the  effect  of  the 
various  magnetic  and  model  parameters  on  them.  Some  elementary 
characteristics  are  also  mentioned.  Conclusions  are  derived 
mathematically  in  Chapters  II  and  IV. 


1,  Variable  inclination  I.  (Figure  3-2). 

The  major  features  of  the  magnetic  anomaly  picture 
due  to  induced  magnetization  in  the  northern  hemisphere  con¬ 
sist  of  a  maximum  on  the  south,  and  a  minimum  on  the  north 
side  of  a  structure.  For  a  symmetrical  model,  if  the  strike 
s  b  90°,  (the  strike  s  is  defined  as  the  angle  between  the 
geological  strike  of  the  two-dimensional  body  and  the  magnetic 
north)  in  the  case  of  I  =  45  °»  the  size  of  the  maximum  and  the 
northern  minimum  are  equal  in  amplitude,  and  the  maximum  and 
minimum  are  symmetrical  about  the  center  of  the  model.  For 
smaller  inclinations  the  northern  minimum  is  the  dominating 
feature.  Figure  3-2  also  shows  that  there  is  a  symmetry  about 
I  s  45-J  in  other  words  for  equal  inclination  differences  from 
I  =  45°  the  amplitudes  of  the  maximum  and  minimum  are  equal, 
and  their  positions  with  the  respect  of  the  center  of  the  body 
are  equidistant. 

The  above  statements  applied  to  a  symmetrical  model 
are  true  for  the  negative  inclinations  (southern  hemisphere) 
except  that  the  north  and  south  sides  of  the  picture  are  inter¬ 
changed. 

The  above  mentioned  symmetry  around  I  **  45°  exists  be¬ 
cause  of  the  symmetrical  model  used  in  the  calculation,  and  is 
due  to  the  change  in  the  relative  importance  of  the  horizontal 
and  vertical  components  of  the  magnetic  field  at  I  =  45°* 

2.  Variable  strike  s  (Figures  3-3  to  3-7)* 

Figures  3-3  to  3-7  show  the  effect  of  the  strike  with 
various  inclinations.  It  can  be  seen  that  with  decreasing  strike 
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angle  the  maximum  of  the  total  intensity  curve  becomes  wider,  it 
moves  toward  the  center  of  the  body,  and  its  minimum  toward  the 
north  decreases*  In  the  s  *  0#  case  there  is  no  minimum  (only  a 
maximum  with  all  the  inclinations)  and  the  maximum  is  symmetrical 
about  the  center  of  the  body.  Because  of  the  type  of  the  model, 
at  I  *  O'  and  s  =  0*  there  will  be  no  anomaly  observed.  The  type 
of  the  model  and  the  increasing  importance  of  the  horizontal  com¬ 
ponent  of  the  magnetic  field  toward  low  inclinations  requires 
that  the  strikes  have  a  greater  effect  on  the  low  inclination 
curves. 

3.  Variable  depth  d  (Figures  3-8  to  3-12) 

4*  Variable  width  W  (Figures  3-13  to  3-15) 

5.  Variable  vertical  extent  D  (Figures  3-16  to  3-18) 

These  three  parameters  of  the  model  were  individually 
changed,  keeping  all  the  other  parameters  constant,  in  order  to 
see  how  the  shape  of  the  total  intensity  curve  respond  to  these 
changes  at  various  inclinations. 

The  variable  depth  curves  show  wider  and  smaller  anoma¬ 
lies  with  increasing  depth.  The  variable  width  curves  cannot  be 
analyzed  alone,  because  the  shape  of  the  curve  depends  on  the  W/d 
ratio.  It  can  be  deduced  from  these  curves  that  when  the  W/d  ?£ 
1,  the  horizontal  extension  (A)  of  the  anomaly  curve  (the  hori¬ 
zontal  extension  at  5 0 %  maximum  value,  or  5 0 i»  minimum  value  at 
low  inclinations)  is  quite  insensitive  to  the  width  of  the  body. 
When  W/d  £  5,  the  change  of  A  almost  entirely  depends  on  the 
width,  and  A  approximately  equals  W. 
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The  flattening  of  the  maximum  peak,  or  at  low  incli¬ 
nation  of  the  minimum  peak,  is  due  to  the  effect  of  the  lower 
edge  of  the  model;  that  is,  the  effect  of  the  vertical  extent 
of  the  model  D«  On  the  curves  near  I  **  45*  there  is  a  shallow¬ 
ing  between  the  maximum  and  minimum  due  to  this  effect  (see  Pig. 
3-17)  which  can  clearly  be  seen  at  smaller  values  of  D.  This 
flattening  effect  is  evident  when  D  £  W.  When  the  width  is  much 
greater  than  the  vertical  extent,  the  maximums  or  minimums  of  the 
anomaly  curve  will  appear  at  the  edges  of  the  body,  while  over 
the  center  the  curve  will  be  shallow,  and  nearly  horizontal. 

This  effect  can  be  seen  also  on  the  variable  width  curves  at  W  = 

15  km,  when  the  width  of  the  model  is  equal  to  its  vertical  extent. 

It  is  evident  that  there  is  an  altitude  of  observation 
above  which  the  effect  of  the  lower  edge  cannot  be  observed. 

6.  Slanted  bodies  (Figures  3-19  to  3-22) 

The  present  calculations  are  over  four  different  models. 
The  tops  are  horizontal  and  symmetrical  about  the  center  of  the 
coordinate  system,  the  sides  are  30*  and  60®  from  the  vertical 
toward  the  north,  and  on  the  other  two  models,  toward  the  south. 

The  lower  edges  are  horizontal  also,  and  they  are  at  15  km  depth. 
The  magnetic  total  intensity  has  been  calculated  for  three  incli¬ 
nations  over  these  models,  in  order  to  show  the  change  on  these 
curves  due  to  the  asymmetry  of  the  models. 

The  general  picture  of  the  magnetic  total  intensity  over 


these  asymmetrical  bodies  is  the  same  as  for  symmetrical  bodies: 
there  is  a  maximum  to  the  south,  and  a  minimum  to  the  north  of  the 
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body  (northern  hemisphere),  but  now  this  maximum  or  minimum  does 
not  depend  entirely  on  the  inclination,  but  depends  upon  the 
shape  of  the  body  also.  The  various  symmetries  in  the  total  in¬ 
tensity  picture,  which  were  mentioned  at  the  group  1  calculations, 
are  only  valid  for  symmetrical  models.  The  total  intensity 
picture  is  different  at  all  inclinations  for  an  asymmetrical  body. 

Because  the  mass  is  more  distributed,  the  amplitudes 
are  smaller;  over  the  slopes  of  the  model  the  anomaly  becomes 
wider,  and  over  the  sharp  tapering  edge  near  the  observation  line 
(surface)  the  anomaly  is  sharp  and  pronounced.  These  effects  can 
be  seen  by  comparing  the  two  oppositely  slanted  models:  Pig.  3-19* 
Pig.  3-21,  and  Pig.  3-20,  Pig.  3-22.  All  four  figures  show  that 
the  shape  of  the  upper  part  of  the  model  has  a  dominating  effect 
on  the  anomaly  curve,  while  the  deep  extension  of  the  model  only 
causes  the  widening  of  the  maximums  or  mlnimums  over  the  sloping 
side. 

Except  for  the  symmetry  relations,  the  same  conclusions 
given  in  the  previous  sections  can  be  drawn  for  these  slanted 
bodies  at  various  strike  angles  s,  and  the  d,  ff,  and  D  changes 
would  also  have  similar  effect  as  they  had  on  the  symmetrical 
models. 
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Chapter  IV 


A  Method  of  Interpretation  for  Anomalies 
of  Total  Magnetlo  Intensity  caused  by 
Two-Dimensional  Bodies 

E.G.  Zurflueh 


A.  General  Remarks 

There  are  numerous  methods  of  Interpretation  for  mag¬ 
netic  data,  each  one  having  a  certain  range  of  application. 

The  method  of  Egyed  (1948)  which  refers  specifically  to  two- 
dimensional  bodies  depends  on  the  measurement  of  the  hori¬ 
zontal  and  vertical  components  and  can,  therefore,  not  easily 
be  used  for  total  Intensity  data* 

Among  the  various  existing  procedures  for  the  analysis 
of  total  Intensity  anomalies,  the  most  widely  used  are  probably 
the  ones  of  Vacquier  et.al.  (193>1)»  Henderson  and  Zietz  (I9f>8) 
and  different  mathematical  methods  as  described  by  Peters  (1949)* 

The  mathematical  methods  are  complicated  and  can  only 
be  applied  to  the  interpretation  of  accurate  magnetic  maps. 

The  more  graphical  methods  on  the  other  hand  have  limitations 
Imposed  by  the  geometrical  configuration  of  bodies  to  be  ana¬ 
lyzed.  The  method  of  Vacquier  et.al.,  for  instance,  becomes 
Inaccurate  when  applied  to  bodies  with  a  horizontal  extent 
equal  or  smaller  than  its  depth  of  burial.  The  method  of 
Henderson  and  Zietz  has  been  partly  devised  to  overcome  this 
difficulty  and  gives  good  results  for  narrow  bodies,  while  it 
is,  in  turn,  not  applicable  to  wide  bodies.  Both  methods  re¬ 
late  to  prismatic  or  cylindrical  bodies  of  more  or  less  iso¬ 
metric  horizontal  cross  seotlon. 
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Since  many  of  the  important  geological  structures  causing 
magnetic  anomalies  are  much  more  extended  in  one  direction  than 
in  others,  the  computations  of  Chapter  III  are  expanded  to 
furnish  a  quantitative  method  of  analysis  for  two-dimensional 
bodies. 

It  is  believed  that  the  method  can  close  a  certain  gap  in 
the  range  of  graphical  methods  and  it  is  also  considered  desirable 
to  work  out  a  simple  and  uniform  method  which  can  be  used  for 
bodies  of  small  as  well  as  large  width  compared  to  depth  of  burial. 

The  method  has  been  devised  for  simplicity  of  procedure. 

In  ord-er  to  minimize  the  influence  of  adjacent  anomalies  only 
distances  near  the  main  extreme  of  the  anomaly  have  been  used  for 
the  depth  determinations.  As  a  general  rule  it  is  only  feasible 
to  use  differences  between  abscissas  of  certain  points  on  the 
curve  since  all  measures  involving  the  total  Intensity  values 
are  dependent  on  the  susceptibility  of  the  body.  Prom  various 
distances  that  have  been  considered,  the  distance  h  between  points 
of  half  maximum  (or  minimum)  value  has  been  selected  as  the  main 
depth  measure  (see  Pig.  4"2)»  The  computations  have  shown  that 
this  distance  (hereafter  called  half -width)  depends  mainly  on 
depth  and  width  and  is  comparatively  little  influenced  by  other 
factors. 

The  distance  a  has  been  chosen  so  that  the  ratio  b  =  h/a 
gives  a  clear  distinction  between  different  values  of  W/d.  In 
order  to  satisfy  this  condition  two  different  distances  a  had  to 
be  used  for  middle  and  for  low  and  high  magnetic  latitudes.  Por 
inclinations  near  0®  or  ±  90°  there  is  a  choice  between  both 
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values  of  a  and  depending  on  the  circumstances  one  or  the  other 
may  give  better  results. 

For  the  practical  applications  it  is  useful  to  keep  the 
limitations  of  a  method  in  mind.  For  this  purpose  the  basic 
assumptions  that  have  been  made  are  briefly  summarized  he  s: 

1.  two-dimensional  body. 

2..  vertical  sides  (this  assumption  implies  that  the 
method  is,  in  general,  to  be  applied  to  anomalies 
which  are  caused  by  susceptibility  contrast  and  not 
by  topography) . 

3.  vertical  extent  larger  than  depth  of  burial. 

4.  uniform  magnetization. 

5.  no  remanent  magnetization  (if  there  is  a  remanent 
magnetization  parallel  to  the  present  earth’s  field, 
the  method  can  be  used  but  susceptibility  estimates 
will  yield  too  high  values). 

B.  Mathematical  Expressions 

The  derivation  of  a  few  mathematical  expressions  for  the 
anomalies  to  be  analyzed  will  be  helpful  in  evaluating  the 
properties  of  the  computed  curves. 

1.  Formula  for  the  anomaly  curve 

Assuming  induced  magnetization  only,  we  obtain  the  following 
expression  for  the  total  intensity  anomaly  produced  by  a  two- 
dimensional  body  with  vertical  sides  and  rectangular  cross  section 
(see  Fig.  4“!)  b7  u&ing  equations  (2-7)#  (2-8),  (2-9)  and  the 
expressions  on  page  2-10! 


4-4 


T  =  Kx  A  +  KgB 


(4-D 


at  which  , 

K,  ~  X  k  F-  C  X  >^s  -  1) 

1^  a  X  k  F  >\/v\  XX  >-Vv\  ,5 
'i. 

““  A  -  e  t  -  e,  -  a,'  +  */ 

B  ^  -  «T  £ 


=  depth  of  burial 
=  width 

=  vertical  extent 
=  field  point  (on 
the  x-axis) 

=  origin  (note  that 
the  origin  used  in 
this  chapter  is  located 
on  the  body’s  axis  of 
symetry  and  not  at  the 
field  points  as  in 
chapter  II.) 


section  of  assumed 


body  and  coordinate  system, 
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Using  distances  instead  of  angles  and  setting  cl  ■  d  +  D  for 
simplification  we  obtain  the  following  expressions  for  A  and 


B: 


^ 

A=  a^Kr^ ~  ^zzk  +  tzs 


■L 


B  ~  ^  OOCy,  :  X 

t  < 


tlVt 


i. 

x.  )  J  ,-°’> 


/  »- 


4 


(4-3) 


V*  *  l  j 


//  i  ■ 

r  '  ,  / 


r  -c 


sr  ■  <■  "i ! 

;r^A  +^'Tyj)  (4-4) 


Equation  (4**4)  can  also  be  written  as: 


1  -  "Ol'-JL'  l ‘  :JV 

(4-4*) 


T* 


Equation  (4—1 )  gives  the  resolution  of  the  total  intensity  anomaly 
into  an  even  and  an  odd  term  as  will  be  shown  here: 

The  function  A(x)  represents  the  even  term  since  replacing 
x  by  -x  in  equation  (4-3)  does  not  change  the  expression.  B(x) 
can  be  written  as  B(x)  =  log  f(x),  where  f(x)  is  the  square  root 
of  the  fraction  in  equation  (4-4*) •  B  is  odd  **  must  satisfy 
the  condition  log  f(x)  =  -log  f(-x),  and  therefore,  f(x)  =  l/f(-x). 
That  this  condition  is  fulfilled  can  easily  be  seen  from  equation 

(4-4*) • 


2.  Extremes  positions 

In  order  to  evaluate  the  maxima  and  minima  positions  we 
take  the  derivative  of  formula  (4-1) 

T»  =  K^A'  +  K2B» 
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Taking  the  derivatives  of  A  and  B  from  equations  (4-3) 
and  (4-4)  w®  obtain: 


3%(v+^f  J\(**!j* 


+(*-&*■ 


1 1 


(4-5) 


*  A. 


jV-:/  Jv<--;V  jW-x/(4-6) 


/ 


In  a  more  concise  notation  the  derivative  of  the  total  intensity 
anomaly  can  then  be  written  as: 

:? 1  .  •-  —i 

(4-7) 


r/  ‘m): 

1  _*-)  - 


where  £  *  •»  I  ( D ,  D  x  D - I>  4  Dj  )  1- "<  ( U,D  ^  u+ -  D.D^  ) 

end  ?'  -  DDjD4)  r  («•- x)( UxbJ 


Expression  (4-7)  is  of  the  form  l  {.<) 


?S) 

r  oc) 


where  p(®)(x)  and  P^)(x)  are  polynomials  of  the  eighth  and  fifth 
degree  respectively. 

Setting  T*  (x)  =  0,  (x)  contributes  the  two  solutions  x  =  ±  w> 

(where  the  anomaly  curve  approaches  the  x  -  axis  asymptotically), 
and  in  the  following  we  will  only  consider  the  equation: 


p(5))x)  =  K]E'  +  K2P«  =  0 


(4-8) 
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Substituting  for  D-j_  through  D|^  we  obtain  the  following  ex¬ 
pressions  for  E»  and  P': 


2  WD  xs  -v  WD(4-  Joi  *  W*-)  y3  (4-9) 

+  W\2JJ(P-S)  -+3^1D-^w4D]x 


FVx')®  3  )(J  +  ‘^"~T  )* 


(4-io) 


For  the  general  case  it  is  not  possible  to  obtain  an  exact  so¬ 
lution  for  equation  (4-8) •  However,  since  E*  is  an  odd  function 
of  the  fifth  degree  and  F'  an  even  function  of  the  fourth  degree, 
we  can  solve  the  problem  for  cases  where  either  or  K£  is  equal 

to  zero. 

a.  K£  *  0. 

According  to  equations  (4-2)  this  condition  is  satisfied 
for  I  e  0*  or  ±  90°  and  for  s  *  0°.  Solving  the  equation  E»  (x)  =  0 
we  get  the  following  extremes  positions: 

x, - O  *w5  =  i /ji  +$l± 

(4-iD 

Id  •  K]_  c 

This  condition  can  also  be  written  as  tan  I  =  ±  sin  s 
(from  (4-2).  The  relation  can  only  be  satisfied  for  +  45^I?-4£0* 

I  =  45®,  a  =  90°  and  1  =  0°,  s  =  0®. 


Two  extreme  values  are: 
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For  the  latter  case  both  Kj  and  £3  are  «qual  to  zero  and  no  a- 
nomaly  will  be  observed. 

From  the  equation  F«  (x)  *=  0  we  obtain  the  extremes  positions: 


for  W^d,  D  the  expression  reduces  to:  x  ■»  ±  W/2  (4-12a) 

3.  Summary 

Formulas  have  been  derived  for  anomaly  curves  over  two- 
dimensional,  rectangular  and  vertical  sided  bodies  and  for  the 
positions  of  the  extremes  of  the  curves  for  special  cases  of 
Inclination  and  strike.  It  is  not  possible  to  obtain  exact  so¬ 
lutions  for  some  other  quantities  like  the  half -width  that  are 
used  In  the  method  of  analysis  presented  here. 

However,  the  formulas  give  information  about  the  general 
properties  of  the  anomaly  curves  and  they  also  allow  conclusions 
about  the  relationship  between  the  different  parameters  Involved. 
Many  of  the  relations  mentioned  here  are  illustrated  in  a  quali¬ 
tative  way  by  the  curves  shown  in  Chapter  III  of  this  report. 

The  following  main  conclusions  can  be  pointed  out: 

The  parameters  k  and  F  are  merely  multiplying  factors  and 
have  no  influence  on  the  shape  of  the  anomaly  curve  (see  equations 
(4-1)  and  (4-2)). 

The  angular  parameters  I  and  s  determine  the  relative  im¬ 
portance  of  the  even  and  odd  terms  of  the  curve.  From  equations 
(4-2)  it  follows  that  for  positive  and  negative  inclinations  of 
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equal  magnitude  the  anomaly  curves  are  symmetric  with  respect  to 
the  z-axis,  if  the  other  parameters  are  held  constant.  This  is 
true  if  the  origin  of  the  coordinate  system  is  taken  over  the 
center  of  the  body,  and  the  profiles  are  drawn  in  the  standard 
way  with  the  northern  end  to  the  right. 

It  can  be  shown  that  for  s  =90°,  constant  body  parameters 
and  for  inclinations  that  differ  from  45°  by  equal  and  opposite 
amounts,  there  is  a  central  symmetry  with  respect  to  the  origin. 
Mathematically  this  can  be  formulated; 

I  «  45°  ±  d  T(x,  ol)  ■  -T(-x,  -cl)  for  s  =  90° 

The  two  symmetry  relations  (see  also  page  2-9)  reduce  the 
amount  of  computation  necessary  for  the  method  of  interpretation. 
Since  we  can  account  for  different  strikes  by  applying  a  cor¬ 
rection  to  the  values  computed  for  s  =  90*  (see  page  4“43)»  we 
can  restrict  ourselves  to  computing  curves  for  Inclinations  be¬ 
tween  0#  and  +  45  °» 

The  first  solution  of  equation  (4-12)  gives  the  distance 
of  the  point  with  maximum  value  from  the  origin  which  coincides 
with  a  zero  position  of  the  anomaly  curve  for  the  cases  where 
Kx  -  0.  This  is  the  distance  a  which  is  used  for  the  interpre¬ 
tation  at  certain  Inclinations  (see  Pig.  4“2).  The  formula 
shows  that  for  great  width  It  is  approximately  equal  to  half  the 
body’s  width  and  that,  therefore,  the  main  extremes  are  located 
over  the  edges  of  the  body. 

With  respect  to  the  body  parameters  d,  W  and  D  it  can  be 
said  that,  as  test  computations  have  shown,  ohanges  In  D  in 


general  affect  mainly  the  amplitude  of  the  anomaly,  while  they 
have  less  influence  on  the  location  of  the  extremes  and  of  the 
zero  positions.  Furthermore,  for  physical  reasons  it  is  clear 
that  if  D  is  greater  than  a  certain  value,  say  greater  than  d 
or  W,  an  increase  in  D  will  only  slightly  affect  the  anomaly- 
curve  because  the  added  magnetized  mass  is  too  distant  from  the 
plane  of  measurement. 

C.  Description  of  Method 

The  quantitative  analysis  of  total  intensity  anomalies  de¬ 
scribed  below  is  based  on  theoretical  curves  which  have  been  com¬ 
puted  for  rectangular  and  vertical  si^ed  two-dimensional  bodies 
(see  Fig.  4“l)»  Sets  of  curves  have  been  obtained  for  fixed  ratios 
of  width  to  depth  by  varying  deoth  and  width  proportionally  (Figs. 
4-11  to  4~26).  In  addition  the  curves  for  variable  strike  of 
Chapter  III  (Figs.  3-3  to  3-7)  are  used.  The  analysis  can  be 
carried  out  in  the  following  steps: 

1)  The  inclination  I  of  the  region  is  determined  by  some 
means  (isogonic  maps,  tables  or  measurements). 

2)  The  approximate  strike  s  of  the  body  causing  an  a- 
nomaly  is  determined  from  the  magnetic  or  geological 
mao  or  from  other  data  available. 

If  no  information  of  this  kind  is  at  hand  a  rough 
estimate  of  s  can  sometimes  be  obtained  by  comparing 
the  measured  anomaly  with  the  curves  with  variable 
strike  in  Chapter  III. 

3)  If  a  magnetic  map  is  to  be  interpreted,  a  profile  is 
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drawn  through  the  center  of  an  anomaly  at  right  angles 
to  its  strike. 

If  only  profiles  have  been  measured,  the  profile 
has  to  be  projected  on  the  direction  perpendicular  to 
the  strike.  For  the  practical  application  it  is  suf¬ 
ficient  to  multiply  the  distances  h  and  a  (see  below) 
by  the  sine  of  the  angle  between  profile  and  strike. 

4)  The  profile  then  is  compared  with  the  computed  curves 
and  a  curve  that  most  closely  resembles  the  actual  one 
is  chosen,  while  keeping  in  mind  the  changes  in  shape 
produced  by  the  strike. 

Since  often  the  exact  regional  field  is  not  known, 
the  main  purpose  of  this  comparison  is  to  determine  a 
zero  line  for  the  anomaly. 

5)  From  the  measured  anomaly  profile  the  dist-nces  h  and 
a  and  the  ratio  b  are  determined  from  the  main  extreme 
of  the  curve  (see  Fig,  4-2).  It  is  essential  that  the 
part  of  the  anomaly  which  lies  between  the  main  maximum 
and  main  minimum  be  used  for  the  determination  of  a. 

h  *  x’  -x  .  (half -width) 
0.5  0.5 

a  =  xo.5  ”xo  or  Vxo 

according  to  the  inclination 
(see  nomograms  A), 
b  =  h/a 

Fig.  4-2:  Total  intensity  anomaly  and  explanation 


of  h,  a  and  b 
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6)  If  the  strike  s  differs  substantially  from  90®  a  cor¬ 
rection  has  to  be  applied  to  h  and  b  by  using  tables 
4-1  to  4-4  (page  4-55). 

7)  With  the  corrected  values  of  h  and  b  enter  nomogram  A 
for  the  nearest  inclination.  (Pigs.  4“2®  to  4-37). 

Read  a  ratio  of  W/d  for  the  body. 

8)  Nomogram  A  also  gives  an  estimate  of  the  depth  of  burial 
d.  A  more  accurate  depth  estimate  can  be  obtained  by 
using  nomogram  B  with  h  and  the  previously  determined 
value  of  W/d.  The  width  W  of  the  body  can  then  be  cal¬ 
culated  from  d  and  W/d. 

9)  For  inclinations  near  the  middle  of  the  l£®  intervals, 
of  the  computed  curves  an  Interpolation  is  necessary. 

For  this  purpose  a  W/d  -  ratio  is  determined  from  nomogram 
A  for  the  two  adjacent  inclinations.  Then  the  average  of 
the  two  W/d  -  ratios  Is  taken  and  a  depth  d  determined 
from  nomogram  B  for  both  inclinations.  The  average  of 
these  two  values  gives  the  depth. 

10)  The  susceptibility  k  of  the  rock  body  can  be  calculated 
from  the  maximum  value  of  the  anomaly.  The  model  a- 
nomalies  have  been  computed  for  k  =  10~3  emu  and  a  total 
field  strength  of  £0,000  gammas.  After  having  corrected 
for  the  total  Intensity  of  the  earth’s  field  in  the  area, 
we  obtain  the  susceptibility  in  electromagnetic  CGS-units 
by  the  following  relationship: 
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k  =  10"3  T  max/T'max 

where  T  max  =  maximum  value  of  measured  anomaly 
T'max  ■  maximum  value  of  computed  anomaly 
It  has  to  be  borne  in  mind,  however,  that  the  vertioal 
extent  D  and  the  strike  s  have  a  strong  Influence  on 
the  amplitude  of  the  anomaly,  and  therefore,  affeot  the 
susceptibility  measurements* 

The  Interpretation  procedure  can  be  clarified  by  stating  an  ex¬ 
ample*  The  profile  in  Fig.  4-3  has  been  computed  for  an  incli¬ 
nation  of  26°  and  the  true  body  parameters  are  indicated  in  the 
figure. 

The  half -width  h  taken  from  the  profile  is  6*5>  km,  tne 
distance  a  «  xm-xQ  is  4*2  km  and  the  ratio  b,  therefore,  equals 
1,55*  Correcting  for  the  strike  according  to  table  4-1  add 
4#  to  the  half -width  and  obtain  h  »  6.8  km.  Similarly  the  cor¬ 
rected  value  of  b  is  1*44« 

Entering  nomogram  A  on  page  4-£l  the  adjusted  values 

of  h  and  b,  we  estimate  a  W/d  -  ratio  of  3.3.  with  this  and 
with  h  nomogram  B  on  page  4 -52  gives  us  the  estimated  depth  of 
the  body  which  is  approximately  2  km.  < 

In  order  to  complete  the  interpretation  it  is  necessary 
to  determine  the  location  of  the  body's  center.  This  can  be 
done  in  a  more  or  less  qualitative  fashion  as  follows: 
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According  to  eq.  (4-11)  for  I  *  0*  or  ±  90*  and  for  s 
a  90*  the  anomaly  is  symmetric  with  respect  to  a  vertical 
line  through  the  center  of  the  body  and  the  main  extreme  is 
looated  over  the  center. 

For  the  cases  where  the  anomaly  curve  represents  an 
odd  function  (in  the  standard  coordinate  system,  see  eq.  (4- 
12).  The  maximum  and  minimum  are  located  approximately  over 
the  edges  of  the  body,  while  its  center  is  given  by  the  zero 
position  between  them. 

For  other  inclinations  and  strikes  the  center  of  the 
body  can  be  found  between  the  main  extreme  and  zero  positions 
by  interpolation  according  to  the  computed  curves. 

D.  Results  of  Practical  Applications 
1.  Analyzed  Cases 

Some  anomalies  on  Canadian  aeromagnetic  maps  have  been 
interpreted  with  the  method  described  above.  From  geological 
map3  it  can  be  inferred  that  the  chosen  anomalies  are  caused 
by  rocks  exposed  on  the  surface.  The  depths  determined, 
therefore,  can  be  directly  compared  with  the  flight  altitudes. 
Since  no  susceptibility  values  measured  from  rock  samples 
were  immediately  available,  no  susceptibilities  have  been  cal¬ 
culated  for  the  examples. 

Fig.  4-4  to  4-^0  show  parts  of  the  aeromagnetic  maps 
with  the  analyzed  anomalies.  The  chosen  profiles  are  indicated 


by  the  lines  AB  and  the  distances  h  and  a  are  also  shown  in 
the  figures.  The  zero  level  used  is  given  by  the  point  0 
on  the  profiles.  The  inclination  and  the  strike  for  each 
case  are  indicated  in  parenthesis  in  the  accompanying  text. 

a.  Port  Hawke sbury,  H.S.  (Pig.  4-4) 

The  ratio  of  W/d  was  found  to  be  approximately  7  and 
the  depth  determined  is  290m  (I  «  73®#  a  “  37*) •  The  flight 
altitude  is  about  1000  feet  or  305m  above  ground.  The  geo¬ 
logical  map  shows  an  outcrop  of  Pre-Carbonlferous  granitic 
rocks  surrounded  by  Mlssisslppian  sediments  at  the  location 
of  the  anomaly. 

b.  Deskenatlata  Lake  South,  NW  Territories 

For  the  anomaly  near  the  Taltson  River  (Pig.  4-5)  a 
W/d  of  8  and  a  depth  of  290m  were  obtained  (I  =  81°,  s  * 

22°). 

The  anomaly  A^  in  Pig.  4“6  yielded  a  W/d  ratio  of 
4  and  a  depth  of  225m  (Ik  8l*,  s  *  13°)#  the  profile  Ag  Bg 
a  W/d  of  2.5  and  a  depth  of  320m  (I  =  8l®,  s  **  10°). 

The  flight  altitude  for  all  three  cases  is  1000  feet 
or  305m.  The  average  of  the  three  depths  determined  is  278m 
which  differs  less  than  10$  from  the  average  flight  altitude. 

Geologically  the  area  is  composed  of  Archaean  and  Pro¬ 
terozoic  acid  rocks. 

Since  75®  Is  the  nearest  inclination  and  because  the 
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graphs  for  I  ■  75®  are  more  reliable  for  email  depths  than 
the  ones  for  I  *  90°,  no  interpolation  was  attempted  and  only 
the  graphs  for  I  «  75®  were  used  for  these  cases. 

c.  Mack  Township,  Ontario  (Pig.  4“7) 

With  a  W/d  ratio  of  2  we  estimate  the  depth  to  be  170m 
(I  *  76°,  s  -  90#).  The  average  flight  altitude  is  500  feet 
or  1 53m  above  ground.  The  rocks  of  the  area  are  Huronian 
quartzites  with  quartz  diabase  and  quartz  norite  intrusions. 

d.  Gladstone  Township,  Ontario  (Pig.  4-8) 

The  estimated  depth  is  130m  and  the  W/d  ratio  1  (I  m 
76*,  s  =  55*) •  The  geological  map  indicates  Huronian  sedi¬ 
ments  of  the  Gowganda  Formation  with  intrusives  of  the  same 
type  as  found  in  Mack  Township. 

e.  Bright  Township,  Ontario 

The  profile  in  Pig.  4“9  gives  a  W/d  of  1,3  and  a  depth 
of  120m  (I  *  76°,  s  *•  90®) •  The  results  for  Pig.  4”10  are: 
W/d  «=  l/2  or  smaller  and  d  *  170m  (I  ■  7 6°,  s  =  90®). 

Pre-Huronian  g ranite  gneiss  and  related  plutonic  rocks 
are  the  rock  formations  encountered  near  the  anomalies. 

The  average  value  of  all  depth  estimates  for  the  ex¬ 
amined  anomalies  in  Ontario  (which  are  relatively  close  to¬ 
gether)  is  l48m,  which  comes  very  close  to  the  average  flight 
altitude  of  153m. 
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2.  Conclusions 

The  test  cases  have  shown  that  the  method  gives  reliable 
depth  estimates  for  anomalies  that  are  reasonably  elongated 
in  shape.  For  carefully  chosen  anomalies  the  results  often 
are  accurate  to  within  10$,  although  only  very  small  depths 
of  burial  have  been  considered.  For  greater  depths  it  is  to 
be  expected  that  the  method  will  give  even  better  results, 
because  deviations  from  the  basic  assumptions  (uniform  mag¬ 
netization  etc.)  will  be  less  important  for  those  cases. 

For  this  method,  like  for  others,  it  is  of  advantage 
to  choose  isolated  and  simple  anomalies.  As  a  general  cri¬ 
terion  it  can  be  said  that  good  results  can  be  expected,  if 
the  shape  of  the  anomaly  profile  is  in  its  major  parts  com¬ 
parable  to  the  theoretical  curves. 

For  anomalies  which  are  nearly  isometric  in  the  map 
projection,  the  depth  estimates  obtained  with  this  method 
tend  to  become  too  small,  because  in  this  case  the  bodies 
considered  are  not  two-dimensional.  However,  if  an  anomaly 
of  more  or  less  elliptic  shape  is  about  three  times  as  long 
as  broad,  the  method  can  be  applied  successfully. 

It  was  not  possible  to  check  the  W/d  ratios  determined 
for  the  actual  anomalies  considered  above.  But  measurements 
on  computed  curves  have  shown  that  they  are  somewhat  less  re¬ 
liable  than  the  depth  estimates. 
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FIG.  4-4=  PORT  HAWKESBURY,  N.S. 
MAP  238 G 


FIG.  4-6’  DESKENATLATA  LAKE  SOUTH,  NW  TERRITORIES 


FIG.  4-8=  GLADSTONE  TOWNSHIP,  ONTARIO 
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FIG.  4-9!  BRIGHT  TOWNSHIP,  ONTARIO 


FIG.  4-10=  BRIGHT 
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E.  Curves  and  Graphs 
ls  Computed  curves 

The  bodies  for  which  the  anomaly  curves  have  been  com¬ 
puted  are  of  the  kind  described  in  paragraph  B  and  the  co¬ 
ordinate  system  used  is  essentially  the  same  as  shown  in  Pig. 
4-1.  As  abscissa  the  ratio  of  the  horizontal  distance  x  to 
the  depth  of  burial  d  for  each  curve  has  been  plotted  in 
order  to  facilitate  the  graphical  representation.  The  verti¬ 
cal  extent  D  is  the  same  for  all  bodies,  namely  15>  km.  The 
depth  of  burial  d  in  km  is  indicated  on  each  curve.  In  all 
cases  the  assumed  total  intensity  P  of  the  undisturbed  field 
is  5 0,000  gammas  and  the  susceptibility  k  is  10"3  emu. 

For  width  to  depth  ratios  smaller  than  one  the  curves 

practically  do  not  change  their  character  except  that  their 

amplitude  Is  multiplied  by  a  factor  proportional  to  W/d. 

That  is  the  curves  for  W/d  »  l/2  e.g.  are  approximately  the 

same  as  the  curves  for  W/d  —  1  multiplied  by  the  factor  of 

l/2.  Therefore,  only  curves  for  W/d  ^  1  are  shown  here. 
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FIGURE  4-12 


J  »  4  1*0*  S-90* 


4-28 


FIGURE  4-13 
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FIGURE  4-14 
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FIGURE  4-15 
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FIGURE  4  —  16 
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FIGURE  4-17 
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FIGURE  4-18 
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FIGURE  4-19 
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FIGURE  4-20 
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FIGURE  4-21 
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FIGURE  4-23 
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FIGURE  4-25 
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FIGURE  4-26 


ayn 


■■■■■ 

■■■■■ 

h! 

m 


ummuiisisn  irn 


s 


a 


aaaaaia 

RSSSSSSSSSSSSSSSSI 

■■—■■■■■■■■■■■■■■■■■■la 


■■■■■■■■■■■■■■■■■■■■■■■■I 

Ississ:  ssssssasssssssssssssssa 

■■■■■■■■ 


: : :  : : ;  I ii  i 


assss 

BBSS 


H  Til  t  t 


Mftirn 


jissasi 

■■■■■■SB 

■■■■■■■■ 


I  H-H+H-HI- 


T-T  1  T-H-r  H 

-t :  :  "i  a  “ 

-rrr  -t-  -t-- 


MM—HW 


Eli 


S*90* 


4-42 


4-43 


2,  Nomograms  and  Tables 

The  values  of  the  tables  and  the  nomograms  have  been 
obtained  by  measuring  the  corresponding  distances  on  the  com¬ 
puted  curves. 

a.  Nomograms 

Nomogram  A  for  each  Inclination  shows  two  families  of 
curves.  Both  of  them  give  the  half -width  h  as  a  function  of 
the  ratio  b.  For  the  solid  curves  the  W/d  ratio  is  held 
constant,  while  the  dashed  curves  are  plotted  for  constant 
depths  of  burial. 

On  nomogram  B  the  half -width  is  plotted  as  a  function 
of  depth  for  constant  W/d  ratios. 

b.  Tables 

As  equations  (i|.-l )  and  (4-2)  show,  the  strike  influences 
anomaly  curves  for  different  body  parameters  in  a  different 
way.  Corrections  for  the  strike  for  different  values  of  depth 
and  width  are  given  in  tables  i|.-l  4”4»  ^  4  and  ^  b  are 

the  percentages  that  have  to  be  added  to  the  measured  values 
of  h  and  b  in  order  to  obtain  the  values  which  would  have  been 
measured  for  a  strike  of  90  *• 

For  inclinations  of  less  than  45°  and  strikes  near  0°  the 
positive  part  of  the  anomalies  becomes  larger  than  the  negative 
one,  and  in  these  cases  it  is  more  convenient  to  measure  the 
distances  h  and  a  from  the  maximum  instead  of  using  the  minimum 
(see  figure  1^-2. 7). 
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FIG.  4-37 
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Table  4-1 

Correction  for  the  Strike:  W/d«=5»  d®2,  D»l5 


8 

4h<*) 

AW) 

AW) 

AW) 

1=15° 

6o° 

+  1.8 

+  2.8 

45° 

+  4.1 

+12.4 

30° 

+10.0 

+26.9 

15° 

+14.4 

+23.4 

+11.4 

+14.4 

0* 

"  4.2 

+15.9 

1*30° 

6o® 

+  3.1 

-  3.3 

45° 

+  4.4 

-10.9 

-  1.1 

-36.1 

30° 

-  1.3 

-33.4 

+  4.4 

-11.7 

13° 

-  3.4 

+  3.5 

0° 

-14.7 

-  7.2 

i=45° 

6o° 

+  3.3 

+13.4 

45° 

0.0 

+21.6 

30° 

-  4.0 

+35.7 

15° 

-13.1 

+33.4 

0° 

-16.5 

+19.7 

i=6o* 

60° 

-  2.7 

+  2.2 

45° 

-  5.8 

+  2.1 

30° 

-  8.5 

+  0.9 

15° 

-14.1 

-  5.7 

0° 

-14.8 

-  7.9 

1=75° 

45° 

-  1.5 

-10.4 

i^° 

-  5.1 

+  2.9 

0° 

-  4.8 

+19.7 
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Table  4-2 

Correction  for  the  Strike:  W/d»5»  d“0.5»  D*l5 


8 

Ah(*) 

Ab(g) 

AM) 

1=15° 

45° 

+2.0 

-  1.5 

15° 

+3.2 

-12.7 

+4.9 

+  2.8 

0* 

-1.2 

+330 

1=30° 

45° 

+4.2 

-4.5 

—6.5 

-44.o 

15° 

-0.8 

+10.7 

0° 

+5*4 

+58.0 

1=45° 

45° 

+6.4 

+34.0 

15* 

-1.6 

+57.0 

0° 

-3.1 

+123 

i=6o# 

45° 

+1.7 

+12.6 

15* 

-3.2 

+18.6 

o# 

-4*3 

+60.5 

1-75° 

45° 

-1.5 

+  2.0 

i5# 

-1.9 

+72.4 

o • 

-2.1 

+172 
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Table  4-3 

Correction  for  the  Strike:  W/d«l,  d«2,  D«l5 


8 

Ah**) 

A5(#) 

A 

I»15® 

•6o# 

-  2.6 

-14.8 

45° 

-  4.8 

-31.0 

30® 

-12.6 

-42.1 

15® 

-31.3 

-72.8 

-30.0 

-71.0 

0° 

+  7.2 

+185 

1=30° 

6o® 

-  4*3 

-  5.8 

45° 

-12.8 

-19.3 

-30.8 

-37.9 

30° 

-28.5 

-37.4 

-14.3 

-26.0 

15® 

+  8.8 

+21.6 

0® 

+20.1 

+190 

i*45° 

6o® 

+  7.4 

+15.0 

45® 

+16.2 

+24.4 

30® 

+30.2 

+55.0 

15® 

+44* 0 

+114 

0® 

+52.0 

+306 

l*=6o® 

45° 

+11.5 

+32.2 

15° 

+20.0 

+99.4 

0® 

+21.0 

+184 

I«75® 

45® 

+  1.9 

+24.5 

15® 

+  3.0 

+91.7 

0® 

+  5.o 

+154 
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Table  4 "4 


Correction 

for  the  Strike 

1:  W/d-l, 

d»o.5,  D-15 

a 

&h(*) 

&b(*) 

1-15  # 

45° 

o 

• 

u\ 

1 

-37.0 

i5° 

-35.6 

-81.7 

-29.2 

-75.8 

o# 

+  4*6 

+4oo 

1*30° 

45° 

— 14«  o 

-23.6 

-34.2 

-41.7 

15° 

+10.6 

+19.6 

0° 

+25.0 

+440 

i*45° 

45° 

+19.2 

+23.8 

15° 

+58.0 

+188 

0* 

460.7 

+730 

i«6o° 

45* 

+12.5 

+23.6 

15° 

+22.7 

+160 

0* 

+23.8 

+396 

1-75* 

45° 

+  1.8 

+39.4 

15° 

+  5*6 

+190 

0s 

+  6.5 

+360 

Appendix 


A  Digital  Computer  Program 
for  the  Computation  of  Magnetic  Anomalies 
over  Two-Dimensional  Polygonal  Bodies 

M.  Talwanl  and  J.R.  Heirtzler 


1*  Framing  the  Problem 

Rewrite  eqe  (2-7)  and  (2-8)  for  a  prism  of  Fig.  (2-4)  expressing 
4>  in  terms  of  distance: 


*.  -  **  «  *  -  X. 


«■  *  "  "*i  *  *  -*»/ 


5  <*V\  2T 


2/ 


,  j>  = 


ri 


(*Z,  +  X\  ) 


<^>S  4>  =  —  *'*  =  *,«. 


*Z/  +  *\ 


K  +  *1 


su 


t  - 


*i/  *-  */\  -ijt.  +  *Z> 


A-l 


A-2 


This  yields 


Vs  Z  (M,  Q'  -  M^P') 


(A-l) 


#  =  Z  (M,  P*  +  H*Q') 


with 


(A-2) 


f  ~  -  —  (e,  -  j  -  ■*  n  */«.  /  /  „  v 

‘  '  n— — -  ^-rf'V.r-.j 


*1,  «-*! 


(A-3) 


(A-4) 


By  refering  back  to  page  2-7  ,  it  ia  seen  that  integrations  leading 
to  this  result  were  carried  out  in  a  counter-clockwise  direction, 
that  is  ©  moving  c.c.w.  during  the  integration.  For  a  computer 
program  there  is  an  advantage  in  using  a  different  labeling  of 
body  points  (see  Fig.  A-l): 


1 


Tho  previous  calculations  can  be  utilized  provided  we  replace 
ri  by  *1  ,  n.  by  r,  ,  6,  by  6%  ,  and  ©,  by  ©,  .  These  changes 
give  the  following  equations: 


*  Q  -  Mt  P) 

(A-5) 

-  2  (Hm  P  +  M9i Q  ) 

(A-6) 

P  = 


(*,-  *>,)  * 


*/ 


f  X 
T  *4 ft. 


*1/  ^/L 


W*. 


(A-7) 


0  =  *il.  * <iL_  (s>t-0J  -  (A— 8) 

A/t  ^  #/il 

For  a  given  field  point  (point  of  observation  or  origin  location) 
the  computer  computes  eqs  (A-5),  (A-6),  and  the  total  intensity 
anomaly  (eq  2-9)  for  each  face  of  the  body  using  eqs  (A-7)  and 
(A-8)  as  internal  subroutines.  It  sums  the  anomalies  for  each 
consecutive  face  and  then  proceeds  to  do  the  same  for  other  bodies 
present.  After  the  various  sums  are  printed  for  a  given  field 
point  it  goes  to  the  next  field  point,  etc. 

The  program  has  been  written  specifically  for  bodies 
with  induced  magnetization  only,  but  can  also.be  used  without 
modification  for  bodies  having  remanent  magnetization  or  bodies 
with  both  remanent  and  induced  magnetization. 


For  computing  convenience  we  can  write  (A-5  and  (A-6) 

above  as: 

V  -  2.kF  [(cot  I'XsU*  s'J  Q  - 

H  rr  z  -k  f[o~o i  s‘)r  -h  isu~  z'scQ 

For  bodies  with  induced  magnetization  only,  k  represents  the 
susceptibility  and  F  the  magnitude  of  the  total  undisturbed  field. 
I'  =1#  s’  c  s  as  in  Fig.  2-7.  For  bodies  with  remanent  or  mixed 
magnetization  kF  has,  of  course,  no  significance;  it  is  merely 
retained  for  convenience  so  that  one  can  use  the  same  program  as 
for  bodies  with  induced  magnetization.  In  practice  kF  is  made 
to  represent  the  magnitude  of  Mpem  for  bodies  with  remanent  mag¬ 
netization  by  giving  k  the  dummy  value  of  1  and  F  the  value  for 
Mrem«  Also  I'  «  a,  s  *  b  as  mentioned  on  p.  2-11.  Similarly  for 
a  body  with  mixed  magnetization  k  is  again  given  the  dummy  value 
of  1  and  F  the  value  of  I*  ®o**,  s*  as  mentioned  on 

P.  2-12. 

We  note  that  the  expression  (2-9)  for  the  evaluation  of 
T,  the  total  anomaly,  remains  unchanged  whether  the  magnetization 
is  Induced,  remanent,  or  mixed. 

2.  Flow  Diagram 

To  assist  in  the  understanding  of  the  flow  diagram  and 
the  subsequent  program  the  following  list  of  symbols  is  defined: 


(A-9) 

( A-10 ) 


Table  A-l 


Symbol 

DIP* 

DIFD* 

D* 

DD* 

F* 

C* 

EXX,  ZEE* 

JTOT* 

CONS* 

PO* 

DELX* 

THETA.,  THETB 
OMEGA 

K 

KTOT* 


Program  Symbol! 


Meaning 

Inclination  or  Dip,  I,  in  degrees 

Modified  Inclination  or  Dip,  I* 
in  degrees 

Strike,  s,  in  degrees 

Modified  Strike,  s',  in  degrees 

Total  Magnetic  Field  intensity, 
in  gammas 

Magnetic  susceptibility,  k,  in  e*m.u. 

Coordinates  of  body  comers,  entered 
in  program  in  clockwise  sequence  in¬ 
cluding  first  twice 

Total  number  of  body  points  including 
first  twice 

Constant  value  of  z  of  field  points, 
zero  for  shipboard  surveys,  negative 
for  aeromagnetic  surveys  taking  sea 
level  as  t  ■  0 

x  value  of  first  field  point 

Increment  in  x  value  of  field  points 

Angle  to  corner  measured  from  x-axls. 

Angular  span  of  face  from  field 
point, 

Subscript  for  consecutive  field 
points 

Total  number  of  field  points 
Subscript  for  body  number 


LNO* 


LNOT* 


Total  number  of  bodies 


H 

V 

PSUM 
Q3  UM 
HSUM, 

HASUM, 

TSUM 

TASUM 


Horizontal  anomaly  for  single  face 
of  body  at  field  point,  in  gammas 

Vertical  anomaly  for  single  face  of 
body  at  field  point,  in  gammas 

The  sum  of  P*s  for  the  various 
faces  of  the  body 

The  sum  of  Q's  for  the  various 
faces  of  the  body 

VSUM  Horizontal  and  Vertical  anomalies, 

in  gammas,  at  field  point  for 

given  body 

VASUM  Horizontal  and  vertical  anomalies, 

in  gammas,  at  field  point  for  all 

bodies 

Anomaly  in  total  intensity,  in 
gammas,  at  field  point  for  given  body 

Anomaly  in  total  Intensity,  in  gammas, 
at  field  point  for  all  bodies 


*Data  which  must  be  entered  into  computer  from  typewriter, 
D  and  DIP  are  used  in  the  calculation  of  V  and  H, 

DD  and  DIFD  are  used  to  oaloulate  TSUM  and  TASUM. 
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I  CALC  P,Q  !■ 


>3jMSN3IpW8| 


ENTER  DD,  DIPD, 
KTOT,  LNOT.  CONS 

~i  - 


400 


ENTER  FO,  DELX 
CALC  PX(K).PZ(K) 


ENTER  LNO,C,F, 
JTOT,D,DIP, 

Ey(j),TTw>/T\ 


105 


PSUM 

Q3UM 


CALC  X1,Z1, 
THETA,  R1 
~1l  - 


CALC  12,  Z2, 
THETB,  R2 


no 


,-T 

-I  Is  Z1 


Z2  ?  L  y®8 


p  =  0 

a  t  ° 


r-i - — 1 

1  SENSESWITCH  3  h 


off 


on 


_ 

DETAILED 

PRINT-OUT 

± _ — - t  . 

PSUM  =  PSUM  +  P 
Q3UM  «  Q3 UM  +  Q 


ADVANCE 
SUBSCRIPTS  FOR 
BODY  CORNERS 


I  LAST  FACE  ?|  ,  B9 

jyea 


CALC  VSUM(k) 

hsum(k) ,tsum(k) 

i  — 


PRINT  LNO,C,F, 
PjPIP  — 


L _  , 

r  36 

T 

VASUM 

HASUM 


)  =  VASUM  (K 
1  =  HASUM (K 


+  VSUM(K 
♦  HSUM(K 


LAST  FIELD  POINT  TO  36 


yes 


LAST  BODY  ?  1—1X0 — *  00  TO  105 


yes 

I  CALC  TASUM  FOR  EACHp 


GO  TO  400 


PRINT  K,FX(K),HSUM(K), 
VSUM(K) ,TSUM(K) _ 
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5  FORTRAN  PROGRAM  FOR  IBM  1620 

DIMENSION  FX(k7) ,FZ(47),VASUM(k7) ,HASUM(k7 
DIMENSION  VSUM (47 ) , HSUM( 47 ) ,EXX ( 30 ) , ZEE ( 30 
400  ACCEPT,DD,DIPD,KTOT.LNOT,CONS 
402  SDD»S IN (.0174533*00) 

CDI  PD**C  OS  (.  0174533*0  IFD } 

SDIFD=SIN ( . 0174533*0 IPD ) 

601  ACCEPT, FO,DELX 

603  DO  604  K“1,KT0T 

H&* 

FX(K)-(FO-DELX)+DELX*RK 

FZ(K)«CONS 

VASUM(K)-0. 

604  hasum(k)**o. 

IO5  ACCEPT,LNO,C,F,  JTOT 
ACCEPT,D,DIP 

410  SD=S IN (.0174533*0) 

CDIP-COS ( ,0174533*DIP) 

SDIPssSIN  ( .  0174533*OIP) 

DO  11  J*l, JTOT 

411  ACCEPT,EXX(  J) ,  ZEE(  J) 

11  CONTINUE 

DO  36  K»1,KT0T 
PSUM®0. 

QSUM-O. 

X1**EXX  (1 )  -FX  (K ) 

Z1*ZEE(1)-FZ(K) 

RSQl=0a**2+Zl**2 
IF(Xl)U0,l40,l80 
110  IF(Z1)120, 130,130 
120  THETA=ATN ( Zl/Xl ) -3 . 1415927 
GO  TO  200 

130  THETA*ATN(Zl/Xl )+3.l4l5927 
go  to  aoo 

140  IF(Z1)150,160,170 
150  THETA*-1. 5707963 
GO  TO  200 
160  THETA=0.0 
GO  TO  200 

170  THETA^l. 5707963 
GO  TO  200 

180  THETA-ATN( Zl/Xl) 

200  J*2: 

201  X2=EXX(J)-FX(K) 

Z2*ZEE( J)-FZ(K) 


RSQa«X2-iH*2+Z2#*a 
IF (X2 )210,240, 280 
210  IF ( Z2>  220 , 230 , 23O  , 

220  THETB*5ATN(Z2/X2) -3. 1415927 
GO  TO  300 

230  THETB^ATN(Z2A2 )+3.l4l5927 
GO  TO  300 
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240  IF (Z2)250, 260, 270 
250  THETB»-i.5707963 
GO  TO  300 
260  THETB=O,0 
GO  TO  300 

270  THETBal.5707963 
GO  TO  300 

280  THETB^ATN(Z2/X2) 

300  IF(Z1-Z2)320,31,320 

31  P-O. 

Q=0. 

GO  TO  32 

320  OMEGAsTHETA-THETB 

IF ( OMEGA ) 3201, 3202, 3202 
32021F(  OMEGA-3 . 1415927  J  330 , 330 , 340 
3201IF ( OMEGA+3 . 1415927 ) 34O, 3 30,330 
330  THETD-OMEGA 
GO  TO  370 

340  IF (OMEGA) 350,360,360 
350  THETD“OMEGA +0.2831853 
GO  TO  370 

360  thetd*omega-6. 2831853 

370  X12=X1-0C2 
Z21-Z2-Z1 
XS  QrfC12**2 
ZSQr&21*#Z 
XZ=Z21*X1 2 

GL»0,5*L0G(RSQ2/RSQ1) 

P»( (ZSQ/(XSQ+ZSQ) )*THETD)+( (XZ/(XSQ+ZSQ) )*GL) 

Q= ( THETIS  (XZ/ (XSQ+ZSQ) ) )-(GL*(ZSQ/(XSQ+ZSQ) ) ) 

32  IF ( SENSESWITOH  3)33,34 

33  H=»2.*C*Ftt(  ( CDIPftSDttP)  +  (SDIP*Q) ) 

V=s2.*C*F*(  ( CDIP&SDttQ)  - ( SDIPfrP ) ) 

PRINT.Xl, X2, Zl, Z2,XSQ, ZSQ, XZ 
PRINT, RSQ1, RSQ2, THETA , THETB 
PRINT, K, J, THETD, GL, P, Q, H, V 

34  PSUM=PSUM+P 
QSUM-QSUM+Q 
X1=0C2 
Z1-Z2 

RSQ1=RSQ2  < 

THETA—THETB 

J*J+1 

JRasJ-1 

IF(JR-JT0T)201,202,202 

202  HSUM(K)«2.*C*F#(  (CDIP*SD#PSTJM)+(SDIP#QSUM) ) 

VSUM(K)»2.*C*F#(  (CDIP*SD*QSUM)-(SDIP*-PSUM) ) 

HASUM(K)*HSUM(K)+HASUM(K) 

VASUM(K)  -VSUM(K )  +VASTJM  (K ) 

36  CONTINUE 

PRINT, LNO,C,F,D, DIP 
DO  37  K=1,KT0T 
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TSUM*  (HSUM(K)*CDIPD»SDD)+ (VSUM(K)*SDIPD) 

37  PRINT,K,FX(K)  ,HSUM(K)  ,VSUM(K)  ,TSUM 
IP (LNO-LNOT) 105, 38, 105 
38  K-l 

381  TASUM«(HASUM(K)*CDIPDttSDD}  +  (VASllM(K)#SDIPD) 
39  PRINT  ,  K ,  FX  (K ) ,  HASUM  (K ) , V ASUM(  K ) ,  TASIJM 
K=*+l 
KR*®-1 

IP (KR-KTOT ) 381, 391*  391 
39I  GO  TO  400 
END 
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3.  A  Simple  iibcample 

Suppose  we  wish  to  determine  the  induced  anomalies  of 
three  field  points  for  an  assumed  body  shape  as  given  in  Pig,  A-2. 


Suppose  the  axis  of  the  body  makes  an  angle  of  b0#  with  magnetic 
north,  the  inclination  is  10®,  the  ambient  magnetic  field  strength 
is  ^0,000  gammas,  the  susceptibility  is  0,001  emu. 

We  have  JTOT  =  4,  LNO  »  1,  LNOT  ■  1,  CONS  =  5.,  KTOT  =  3# 

PO  =  5.,  DELX  =  3. 

Information  in  the  following  form  must  be  supplied  to 
the  machine  by  typewriter 

60.  10.  3  15. 

5.  3. 

1  .001  50000.  4 

60.  10. 

8.  8. 

8.  6. 

11.  6. 

8.  8. 
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With  senseswitch  3  off  numerical  results  will  be  printed  as 
follows : 


(LNO) 

(c) 

(P) 

(D) 

(DIP) 

1 

1.0000000E-03 

50000.000 

60.000000 

10.000000 

(K) 

PX(K) 

HSUM(K) 

VSUM(K) 

TSTJM(K) 

1 

5.0000000 

11.142248 

8.2548*30 

10.93631* 

2 

8.0000000 

-25.503102 

53.009137 

-12. *4585* 

3 

11.000000 

-12.179321 

-43.110189 

-17.873370 
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4*  Console  Operation  for  Lamont  IBM  1620 


Switches: 


Clear  Memory: 


Load  Tape: 


Load  Data: 


Senowiteh  3  ON  for  detailed  print  out  (used 
for  trouble  shooting  only) 

Other  senseswitches  OFF 
Program  switches  OFF 
Other  switches  ON 

Reset,  Insert,  160001000000,  Release,  Start, 
(after  a  few  seconds)  Instant  Stop. 

Mount  object  tape,  set  tape  console  switch 
to  Reel.  Reset,  Insert,  360000000300,  Re¬ 
lease,  Start.  If  machine  falls  to  read,  re¬ 
peat  previous  instructions.  If  functioning 
properly  machine  will  pause  2/3  way  through 
reading  to  think  but  will  proceed  to  end  of 
tape  when  it  will  type  "Load  Data". 

Set  typewriter  tabs  at  approximately 


16  28  45  $Q  70. 

Reset,  Insert,  4907500,  Release,  Start 
Numbers  must  be  typed  in  upper  case. 

Minus  signs  and  decimal  points  in  lower  case. 
Zero  in  either  case. 

A  space  must  be  typed  between  successive  data 
within  a  line  of  data. 

A  record  mark  after  the  last  datum  of  a  line 
is  optional. 


Typing  Error: 


Stop  Output: 


Check  Stops: 


A  decimal  point  is  required  after  floating¬ 
point  data. 

A  decimal  point  is  optional  after  fixed- 
point  data. 

If  incorrect  data  has  been  entered  in  machine 
Reset,  Insert,  4907500,  Release,  Start,  re¬ 
enter  all  data. 

If  incorrect  data  has  been  typed  but  not 
entered  into  machine  Senswitch  4  ON,  Release, 
Start,  Senswitch  4  OFF,  repeat  line  of  data. 
Depress  Stop.  Neatness  can  be  achieved  by 
depressing  Stop  while  carriage  of  typewriter 
is  in  return  motion  following  typing  of  last 
line  of  desired  data.  If  machine  is  not 
stopped  it  will  automatically  return  to  the 
stage  where  it  is  again  ready  to  accept  all 
data. 

If  a  check  stop  occurs  during  data  entry  it 
may  be  due  to  a  parity  error  in  that  part  of 
the  memory  related  to  data  storage:  Reset, 
Insert,  4907500,  Release,  Start,  re-enter 
all  data. 

If  a  check  stop  occurs  at  any  other  time  or 
if  the  above  procedure  fails  clear  memory, 
re-load  tape,  re-enter  all  data. 
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